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ABSTRACT 


An  analytical  treatment  of  the  performance  of  a  gas  gun  has  been  carried  out  by  means  of 
a  modified  “Lagrange  approximation.”  The  actual  local  cross-sectional  area  of  the  high- 
pressure  vessel  of  the  gas  gun,  larger  than  the  area  of  the  base  of  the  gun’s  barrel,  has 
been  taken  into  account  and  its  effect  on  the  performance  determined. 
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SECTION  I 
INTRODUCTION 


In  a  special  materials  test  where  it  was  necessary  to  accelerate  a  pellet  to  moderately 
high  velocities,  well  below  the  speed  of  sound,  a  gas  gun  utilizing  compressed  air  was  used. 
The  projectile  in  this  case  was  a  piston  into  Ihe  forward  face  of  which  the  pellet  was  implanted. 
The  target  was  placed  close  to  the  end  of  the  barrel  of  the  gun. 

The  piston  velocity  was  measured  at  the  target.  On  some  occasions  the  time  between  the 
start  of  the  piston  and  its  arrival  at  the  target  was  also  measured.  To  obtain  a  check  on 
velocities  as  well  as  to  determine  the  necessary  design/ characteristics  for  various  piston 
velocities,  it  became  desirable  to  carry  out  an  analysis  of  the  performance  of  gas  guns, 
even  if  this  could  be  achieved  only  approximately. 

In  the  following  analysis,  friction  between  the  piston  and  the  barrel  walls,  leakage  of  the 
high  pressure  air  around  the  piston,  and  pressure  buildup  in  front  of  the  moving  piston  are 
not  considered. 
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section  n 

ANALYSIS 


The  interior  aerodynamics  of  the  gas  gun  in  the  present  case  may  be  treated  as  that  of  the 
“Lagrange  Problem”  of  Interior  Ballistics,  since  the  velocities  involved  stay  well  below  the 
speed  of  sound.  The  following  simplifying  assumptions  are  made:  a)  the  gas  expands  adia- 
batically;  b)  the  flow  is  one -dimensional;  and  c)  the  velocity  of  the  gas  at  any  instant  varies 
in  proportion  to  the  distance  from  the  closed  end  of  the  pressure  chamber  and  inversely 
with  the  local  cross  section  of  the  chamber.  The  gas  velocity  is  zero  at  the  closed  end  of 
the  vessel  and  equal  to  the  piston  velocity  at  the  base  of  the  piston. 

The  complete  analysis  is  contained  in  Reference  1.  Here,  only  those  parts  of  it  necessary 
for  the  application  of  the  charts  and/or  tables  (see  Appendices)  are  repeated. 

For  a  gun  having  a  chamber  design  such  as  d  in  Figure  1,  the  local  gas  velocity  at 
position  £  is 

£  A  m 

-fVMO  *({)*« 


The  kinetic  energy  of  the  air  moving  in  the  vessel  and  in  the  barrel  behind  the  piston  which 
may  have  arrived  at  a  position  x  after  the  start  is 


,L„+  x 


ir  =  f°  -f  *<«)>'£  "f* 


Po  V„ 


2  V0+Ax  Uq 


A2  £  f 

L„+x)2  V*  J0 


Lo  + 


A(£) 


(2) 


The  kinetic  energy  of  the  moving  piston  of  mass 


W 


g 


is 


.  w.  Ll 

•piston  Q  2 


(3) 


The  external  work  produced  by  the  adiabatic  expansion  of  the  air  initially  in  the  vessel, 
when  the  piston  has  moved  into  position  x  is: 


.xt  l  »lf  v  *  'xh'  (' ‘  )  [ 1 '  Ci  Jt/L,- Pjp )  J 


(4) 


and  is  equal  to  the  kinetic  energy  of  the  moving  piston  and  the  gas.  The  velocity  of  the 
piston  along  its  path  is  then  (for  any  geometry  of  the  pressure  vessel),  after  some  transforma 


tion. 


vx  = 


2  Po  i 

r-\p0 

t'-f) 

[•-( 

i-  po/P  Y"] 

i+  «/L0-/o//3  I  1 

i 

W/4,  I  L0+x 

/>ovo  3  L.  +  « 

c 

1/2 


(5) 


2 


AFFDL-TR-66-159 


The  time  t  needed  for  thepiston  to  travel  a  distance  x  from  the  start  is: 


(6) 


SECTION  III 


NUMERICAL  RESULTS 


The  contents  within  the  brackets  in  the  denominator  of  the  equation  for  vx  depend  entirely 
upon  the  shape  and  size  of  the  high  pressure  vessel.  The  vessel  proper  may  have  a  simple 
geometric  shape  (Figures  la,  b,  and  c),  or  it  may  be  a  composite  of  several  simple  shapes  of 
overall  length,  L0  (Figure  Id).  Its  dimensions  will  be  referred  to  the  radius  r  of  the  barrel. 

For  a  frustrum  (Figure  la),  the  integral  (from  the  denominator  of  Equation  5) 


r  s  f'  ($/Lo)8d(C/Lo) 


Jo  A  ( £ )/  A 


(7) 


assumes  the  form 


r  -  K  . 


(8) 


If  R*  =  R  or  k'  -  k  ,  the  frustum  becomes  a  cylinder,  and 


(9) 


K 
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For  a  rotationally  symmetrical  ellipsoid  (Figure  lb) 

2  t 

r 

0 


7T 

l-T?~ 


lL  rl°  ut  .  (j_\z  — *—= 

Jo  7t(2]£-(Jl)2£)  1  +yi-(r/b)* 


I  +  ^/l-(r/b)*  (r/b)2 


-I 


>/£l 


i+ 


tis. 


-I 


.  when  b/r  = 


(10) 


The  ellipsoid  becomes  a  sphere,  when  a  =  b  =  R.  With  R  =  )3r  the  formula  for  the  sphere  is 
the  same  as  that  for  the  ellipsoid. 

For  a  paraboloid  (Figure  lc) 

T  .  Ell  =  _ ! _  r  Z_LR hi  \y/_R_V _ ! _  _3_1 

V  1  7T2p(l0+€-0  ( R/r)2  - 1  [  l(R/r)2-l  J  Av  r  )  (iyr)2-l  2] 


s£*7  [2("ffT7)^"^T"T]’wh,n  R/r  =  £ 


For  these  three  shapes  the  numerical  values  of  I  have  been  computed  and  plotted  (Figure  2) 
as  a  function  of  the  characteristic  dimensions. 

For  a  composite  vessel  such  as  Figure  Id,  consisting  of  half  an  ellipsoid,  a  cylindrical  part, 
a  frustum,  and  a  short  cylindrical  part  of  barrel  cross  section,  the  volume  is 


V  * 


-€•>-« + *‘f+\  +f*t*l\*  -  v-  «>•<* 


a+^+J,  +  X2 

+  r  J 

a  +  V*. 


(II) 


+  r  1  +a  +i° +4 +  4 ) }  ■ '  r  3  Mt1*  i(v * ) 


-20  (- 


<M+X0)2  X, 


)  +  i^  +  ,  +  X0  +  ^-+X2} 


4 


AFFDL-TR-66-159 


when 

R/r  =  b/r  =  k=/3,  OA  *  fr,  /.A  ■  Xo.  V'=  Xl  .  *2  A  *  X2  1,21 


Then  L  can  be  computed  according  to  the  definition 
e 


V 


o 


(13) 


■°  L  r-tir"10 


The  integral  I  for  the  composite  vessel  becomes 

'MUW?  •  ».».  W  <£  [*•<■■]*  »■' 

+  _^L_]  +(/i  +  X0+X,)t  x2  +(^  +  VXl}  x2  +  -^§—  +(^3 

-  z  (m+v^<a+vx.>  '^vviirTv)! } 


A  numerical  computation  of  the  velocity  »,  of  the  piston  and  its  ‘ravel  time  alon* ^,U.egun 
barrel  has  been  carried  out  for  a  gas  with  the  specific  heat  ratio  y  -  const  1.4.  To  m 
the  results  as  general  and  as  convenient  as  possible,  certain  magnitudes,  which  can  be 
determined  beforehand  from  the  design  values,  have  been  considered  as  parameters. 


ALo 

a  =  “y —  = 
vo 


Lo 

La 


_ length  of  the  vessel _ , 

equivalent  length  of  the  vessel 


V(/q  _  mass  nf  the  piston 
p  v  =  mass  of  compressed  air 


8  =  I  -  31, 


(15) 
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speed  of  sound  in  an  ideal  gas  at  the  initial  pressure  and  density, 

(16) 


311(1  'Z'  =  Pq/@’  versus 


(17) 


rj  =  ~r —  =  piston  travel  ^  the  indepen(jent  variable.  (18) 

'  L0  length  of  vessel 

Using  these  definitions,  the  dimensionless  piston  velocity  and  the  travel  time  may  be  written 
as: 


2(1-'/')  1 

H 

f  l " 

*  Y  1 

r  (x- 1)  1 

{l+aV-'I'J  J 

cu  + 

1  +V 

l  +a  rj  \ 

r,  s  l 

V  (i  +17)3  j 

(19) 


and 


r 


respectively. 


(20) 


v  and  r  ,  versus  rj ,  have  been  computed  for  selected  values  of  the  parameters  a  ,  ou,  S,  and 
\j/,  and  are  given  in  the  tables. 

Because  of  the  assumptions  involved  caution  must  be  observed  in  the  use  of  the  tables.  The 
ratio  v  x  ^  0  o  ls  no1  a  true  Mach  number  although  it  is  the  ratio  of  the  piston  velocity  to 
the  velocity  of  sound  in  the  compressed  gas  in  the  pressure  vessel  at  initial  conditions.  The 
local  speed  of  sound  in  the  air  after  adiabatic  expansion  in  the  barrel  is  always  lower  than 
aQ  and,  therefore,  the  true  local  Mach  number  of  the  air  in  the  immediate  vicinity  of  the 
piston’s  base  is  higher  than  the  ratio  vx  /  o0  . 

The  computation  in  this  report  is  based  upon  the  assumption  that  the  velocities  involved  are 
well  below  the  speed  of  sound.  We  arbitrarily  restrict  the  applicability  of  the  derivation  for 
the  computation  of  gas  guns  such  that  the  local  Mach  number  in  the  gas  is  not  to  exceed  0.5. 
The  maximum  local  Mach  number  is 

■■ft-)  ■&)&)• 


The  local  speed  of  sound  is 


<0.5,  and  v  <  0. 5 


a 

ao 


(21) 
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and  according  to  definition 


°o 


= 


Y9  pv0 


(23) 


with 


-a.  =  (v<?-b  Y 

po  \  v-b  / 


(24) 


(25) 


Vo  +  Ax 


o  _  / v  /vp-b  \y  v  ~  _  / V0+A x  /  Vp  -  m//3  \X _ 

°o'Vv0l  V-b  )  V-b  'V  vo  VVo+Ax-m//3/  V0+Ax-m//3 

/  (l+x/Le)2  /  <-(pp/j9)  \y  _  /(i+QT7)g  / 

-  Vl+(x/Le)  -(/50//5)  Vl+(VLe)Lf>e//3)J  'Vi+aij-^  V 


(26) 


l  +  a  17  -  \f/ 


All  of  the  values  of  v  exceeding 


0.5 


1  1  - 1  Zv 

\l  +  a  17-  \p  } 


(27) 


should  not  be  used.  These  values,  therefore,  have  been  eliminated  from  the  tables  by  the 
expedient  of  printing  zeros  instead.  The  values  of  r  corresponding  to  those  inadmissible 
values  of  v  have  also  been  eliminated  from  the  tables. 
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SECTION  IV 
SUMMARY 


Numerical  values  which  can  be  used  in  the  calculation  of  piston  velocities  and  travel  times 
in  the  barrel  of  a  gas  gun  have  been  calculated.  The  effects  of  different  shapes  of  the  pressure 
chamber  have  been  taken  into  account.  Piston  velocity  and  travel  time  are  given  in  dimension¬ 
less  form  as  functions  of  the  dimensionless  piston  travel  with  certain  design  characteristics 
as  parameters.  These  may  be  computed  beforehand  and  used  for  interpolation  of  the  tables. 

It  is  assumed  that  air  is  the  driving  medium.  However,  the  tables  are  also  valid  for  any  gas 
with  y  =  const  =  1.4. 

In  the  case  of  a  gas  with  a  specific  heat  ratio  other  than  1.4,  the  tabulations  of  the  computer 
programs  are  attached  so  that  the  desired  numerical  values  may  be  computed. 

When  the  pressure  chamber  is  simple  in  shape,  as  illustrated  in  Figure  la  to  lc,  the 
numerical  values  of  I  in  the  parameter  8  =  1-31  can  be  taken  from  the  graphs  in  Figure  2. 


REFERENCE 


1.  W.  F.  Knackstedt  and  J.  L.  Politzer.  “Influence  du  Chambrage  en  Balistique  Interieure” 
(Chambrage  Effect  in  Interior  Ballistics)  (in  French)  Mem.  Artillerie  Franc  (1967) 
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Figure  1.  Representative  Shapes  of  Gas  Gun  Pressure  Chambers 

a.  Frustum 

b.  Rotationally  Symmetrical  Ellipsoid 

c.  Rotationally  Symmetrical  Paraboloid 

d.  Composite 
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Fortran 


APPENDIX  I 
IV  Computer  Program 
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C  FORTRAN  IV  PROGRAM  FOR  COMPUTATION  OF  VELOCITY 


DIMENSION  VFL(ll) 
REAL  MACH 
READ!5*47) GAMMA 

47  FORMAT ( 1F1 2«4 ) 

48  READ ( 5 ♦ 49 ) PS  I 

49  FORMAT! 1F12. 4) 

4  DDELT A  =  .25 


DOMFGA= • 5 

DFTA=.25 

DALPHA=.2 

C  THESE  ARE  STEP  SIZES  FOR  DELTA*OMEGA,ETA*ALPHA 

MACH  =  .5 

C  MACH  IS  LOCAL  MACH  NUMBER 

G  =  GAMMA  -  1. 

F  =  1.  -  PSI 
W  =  2.»U.-PSI)/GAMMA/G 

C  THFSE  ARF  I NTFRMEDI ATE  CALCULATIONS  TO  REDUCE  COMPUTING  TIME 

ALPHA= • 2 
D078  M= 1 ♦ 5 

C  THIS  LOOP  INCREMENTS  ALPHA 

R  =  ALPHA/3. 

C  THIS  IS  AN  INTERMEDIATE  CALCULATION  TO  REDUCE  COMPUTING  TIME 

DFLTA=0. 

0070  L=1 .5 

C  THIS  LOOP  INCREMENTS  DFLTA 


FT  A=  »  25 

WR  I  TF( 6, 40) DFLTA, ALPHA, PS  I 

40  FORMAT! 1H1 ,54X,13HNU  VERSUS  E T A, //27X , 8HDELTA  =  , F4. 2 , 16X , 8HALPHA 

1=  ,F4.2,16X, 6HPS I  =  ,F4.2,//57X  ,5H0MEGA , 

8  //2X,3HETA,6X,3H0.0»8X»3H0.5,8X,3H1 .0»8X, 

23H1 .5,8X,3H2.0,8X,3H2.5,8X,3H3.0,8X,3H3.5,8X,3H4.0,8X,3H4.5,8X,3H5 
3.0) 

DO  38  K= 1 »  20 

C  THIS  LOOP  I NCRFMENT S  ETA 

T  =  1.  +  ETA 
S=  ALOHA*  ETA 
U  =  1.  +  S 
A  =  F+S 

C  THESE  ARE  INTERMEDIATE  CALCULATIONS  TO  REDUCE  COMPUTING  TIME 

OMEGA=0. 


DO  31  1=1,11 

C  THIS  LOOP  INCRFMFNTS  OMEGA 
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VEL(I)  =  (W*(1.-(F/A)**G)/(0MFGA  +  R* T /U* ( 1 . -DFLT A / T**3 . ) ) ) ** . 5 

TFST  =  MACH*(U**2./A*(F/A)**GAMMA)**,5 
I F ( VEL ( I )  -  TEST)  31.31.3 

C  THIS  TESTS  FOR  VELOCITIES  GREATER  THAN  MACH  .5 

3  VEL ( I )  =  0. 

C  THIS  PREVENTS  PRINTOUT  OF  MEANINGLESS  VALUES 

OMrGA=nMrGA  +  OOMrf;A 
WPITF(6*35)FTA,(VFL( T ) ,1  =  1  til ) 

35  FORMAT (1  HO. 1 X. F4*2 *  1 X. 1P1 1  FI  1.3 ) 

38  ETA=ETA+DETA 
70  DELTA=DELTA+DDELTA 
78  ALPHA= AL PH A+D ALPHA 
GO  TO  48 
END 


FIRST  DATA  CARD  IS  FOR  GAMMA 

THE  SUBSEQUENT  DATA  CARDS  ARF  FOR  PST 
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C  FORTRAN  TV  PROGRAM  FOP  COMPUTATION  OF  T I ME 


COMMON  ALPHA*OFLTA,OMEGA,F*G,R*W*GAMMA.PSI*c 
DIMENSION  COR( 21 ) ,A ( 10,21 ) 

REAL  MACH 
EXTERNAL  VEL*SER 

MACH  =  .5 

C  MACH  IS  LOCAL  MACH  NUMBER 

RF AD  ( 5  *19)  GAMMA 
19  FORMAT! 1  FI  2.4) 

13  RFAD(F,18)PST 
18  FORMAT ( F 1 2 • 4 ) 


14  G=  GAMMA  -  1. 

F  =  1.  -  PSI 
W  =  2 • *F/G AMMA /G 

C  THESE  ARE  INTERMEDIATE  STFPS  TO  REDUCE  COMPUTING  TIME 

DALPHA= • 2 
DDFLT  A  =  .23 
DOMFGA= . 3 

C  THFSE  ARE  STFP  STZFS  FOR  ALPHA . DFL T A , A ND  OMFGA 

ALPHA  =  . 2 
D080  L= 1  *  5 

C  THIS  LOOP  INCREMENTS  ALPHA 

R  =  ALPHA/3. 

DELTA=0. 

00  67  K  A  =  1  ,S 

C  THIS  LOOP  IMORFMENTS  DELTA 

WR  I  TE ( 6  *40 ) DELTA ♦ ALPHA , PS  I 

40  FORMAT! 1H1 *60X,15HT IMF  VERSUS  ETA  , //35X , 8HDELTA  =  ,  E4 . 2 . 1 6X , 8HALPH 
1A  =  *F4.2* 16X.6HPSI  =  .F4.2.//65X. 

1 5HOMFGA ,//llX*3HETA»6X,3H0.3*8X, 

23H1.0*8X*3H1.5*8X*3H2»0»8X*3H2.5*8X*3H3.0*8X*3H3.5*8X*3H4.0*8X» 

33H4.5»8X,3H3.0) 

OMEGA= • 5 


C 


DO  51  1  =  1  *  10 

THIS  LOOP  INCREMENTS  OMFGA 


C 


DO  25  K=  2*21 

THIS  LOOP  PERFORMS  INTEGRATION  UP  TO  FTA=5.0 
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K0=K-1 

COR  (  K  )=.25*FLOAT  (  KO  ) 
E  =  COR(K) 


THIS  COMPUTES 

THF 

LOWFR 

LIMIT 

OF 

EACH 

INTEGRAL 

FM=F+.2F 

THIS  COMPUTE 

THF 

UPPER 

LIMTT 

OF 

EACH 

T  NTFGRAL 

H  =  ALPHA  *  F 
Z  =  F  +  H 
Y= 1 • +H 

C  THESE  ARE  INTERMEDIATE  STEPS  TO  REDUCE  COMPUTING  TIME 

TFST=1./VEL(F) 

TACH  =  MACH»Y* (1 . /Z*( F/Z  > **GAMM A ) **. 5 
IF ( TEST-TACH ) 76*76*2 
2  A ( I ,K  )  =  0. 

GO  To  ?  k 

C  THIS  TFFTC  POP  VELOCITIES  LARGFR  THPN  MACH  .5  AND 

C  DPFVFNTS  COMPUTATION  OF  MFAN I NGLFSS  RESULTS 

76  IF  ( K  -  2)8,8,77 

8  A ( I ,K ) =  QUAD! l.E-05 ,l.E-08, 850, SER )+QUAD( 1 . E-02 » 1 . F-05  , 1800, SER 
1 )+QUAD( .25 ,1 .E-02,330»VEL ) 

THIS  PERFORMS  INTEGRATION  BETWFFN  ETA  =  l.F-08  AND  .25 
l.F-08  WAS  SFLFCTED  AS  COMPROMISE  BETWEEN  ACCURACY ,COMPUTAT ION 
TIMP  PFOIII RPMFNTS,  AND  POSSIBLE  OVERFLOW  CONDITION 

GO  TO  28 

7  7  a  (  I  ,K  )  =QtlAD(  FN*F  ,  I  5  ,  VPL  )  +  A  (  I  ,KO) 

C  THIS  INTEGRATES  FROM  .25  TO  5.0  AND  SUMS  UP 

25  CONTINUE 

51  OMEGA  =  OMEGA  +  DOMEGA 

WRITE (6,45) (COR(K) ,<A(I,K),I=1,10),K=2*21) 

45  FORMAT ( 1 HO , 1  OX , 0PF4 . 2 , 1 X , 1 P 1 OE 1 1 . 3 ) 

67  DFLTA=DELTA+DDFLTA 
AO  alpha=al pha+dal°ha 
GOTO  13 
PND 
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FUNCTION  VEL (ETA) 

THIS  IS  THF  RFCIPROCAL  OF  THF  EXPRESSION  FOR  THE  VFLOCITY 

COMMON  ALPHA  , DELTA  »OMFGA  * F  »G » R  * W»GAMMA  »PS I  *  F 

T=1 ,  +  FTA 
S=ALPHA»ETA 
U=1 .+S 
A  =  F  +  S 

THESE  ARE  INTERMEDIATE  STEPS  TO  REDUCE  COMPUTING  TIME 

VEL  =  ( W* ( 1,-{F/A)**G)/ 

1 (OMEGA+R*T/U*< 1 ,-DELTA/T**3. ) ) )#*(-. 5) 

RETURN 

PND 


FUNCTION  SER(FTA) 

THIS  IS  A  SERIES  EXPANSION  OF  VEL  FOR  SMALL  ARGUMFNT 

COMMON  ALPHA, DELTA, OMEGA, F,G,R,W, GAMMA, PSI ,E 

S=ALPHA*ETA 
T= 1 ,  +  ET  A 
U=1 .+S 
T  A  =  S  /  F 

THFSF  ARF  INTERMEDIATE  STFPS  TO  REDUCE  COMPUTING  TIMF 

SFR=(W*G*TA»(1  .-GAMMA/2. »TA*(  l.-(C,AMMA  +  l  .  )  /3.*TA  )  )  / 

1 (OMFGA+R*T/U* ( 1 .-DFLTA/T**3. ) ) )#*<-. S ) 

RETURN 

END 
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FUNCTION  QUAD ( BF,AE*ME* DUM ) 

C  CONTR I RUTFD  RY  JAY  L.  POLITZER 

C  QUAD  PERFORMS  SIMPSONS  INTERGRATION 

BE  =  UPPER  LIMIT 
AE  =  LOWER  LIMIT 
ME  =  NUMBER  OF  INTERVALS 

DUM  IS  DUMMY  FUNCTION  NAME  OF  FUNCTION  TO  BE  INTEGRATED 
REAL  NAME  APPEARS  IN  EXTERNAL  TYPE  STATEMENT ( VEL »SER ) 

COMMON  ALPHA, DELTA, OMEGA, F,G,R,W, GAMMA, PSI ,E 

KE  = ( ME+1 ) /2#2 

THIS  INSURES  PROPER  NUMBFR  OF  INTFRVALS 

DX=  ( BC-AF ) /FLOAT ( Kc ) 

QUAD=DUM ( AE ) +DUM ( BE ) 

EX=AE+DX 
DO  1  1  =  2, KE 

C=4-2*( (10*I/?-I/2*10)/5) 

THIS  PROVIDES  CORRECT  FACTORS,  X.  AND  4  RFSPECTIVELY 

QUAD=C*DUM ( EX ) +QUAD 
1  FX=EX+DX 

OUAD=OUAD*DX/3. 

RETURN 

RND 


DATA 

THE  FIRST  DATA  CARD  IS  FOR  GAMMA 
THE  SUBSEQUENT  DATA  CARDS  ARE  FOR  PSI 


17 


AFFDL-TR-66-159 


APPENDIX  II 
Nu  Versus  Eta  - 

“Piston  Velocity/ Piston  Travel  Distance”  - 
A  Dimensionless  Quantity 
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